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Abstract 

In this study we attempt to quantify the amount of S, Cl and F released by the 1300 km 3 Roza member (~ 14.7 Ma) of 
the Columbia River Basalt Group, which was produced by a moderate-size flood basalt eruption in the mid-Miocene. Our 
results are the first indication of the potential atmospheric S0 2 yield from a flood basalt eruption, and indicate the 
mechanism by which flood basalt eruptions may have seriously affected the environment. Glass inclusions in phenocrysts 
and quenched glass in products from various stages of the eruption were analyzed for concentrations of S, Cl and F and 
major elements. Glass inclusions contain 1965 ± 110 ppm S, 295 ± 65 ppm Cl and 1310 ± 1 10 ppm F. Groundmass glass of 
Roza dike selvages contains considerably lower concentrations: 1110 + 90 ppm S, 245 ± 30 ppm Cl and 1020 ± 25 ppm F. 
Scoria clasts from near vent deposits contain 665 ± 75 ppm S, 175 + 5 ppm Cl and 950 ± 20 ppm F, and the groundmass 
glass of lava selvages contains 520 ± 30 ppm S, 190 ± 30 ppm Cl and 890 ± 55 ppm F. In crystalline lava, the 
concentrations are 195 ppm S, 100 ppm Cl and 830 ppm F. Volatile element concentrations in these samples represent the 
progress of degassing through the eruption and can be used to estimate the potential amount of the volatiles S, Cl and F 
released by the magma into the atmosphere, as well as to evaluate the amount liberated by various phases of the eruption. 
The total amount of volatiles released by the Roza eruption is estimated to have been ~ 12,420 Mt S0 2 , ~ 710 Mt HC1 and 
- 1780 Mt HF. The Roza magma liberated - 9620 Mt S0 2 (77% of the total volatile mass released), ~ 400 Mt HC1 (56%) 
and ~ 1450 Mt HF (81%) at the vents and lofted by the eruption columns to altitudes of 7-13 km. Degassing of the lava is 
estimated to have released an additional ~ 2810 Mt S0 2 , ~ 310 Mt HC1 and ~ 330 Mt HF. 

The Roza eruption is likely to have lasted for ~ 10 years, indicating an annual H 2 S0 4 -mass loading of ~ 1800 Mt. 
Thus, the atmospheric perturbations associated with the Roza eruption may have been of the magnitude predicted for a 
severe “nuclear” or “volcanic” winter, but lasting up to a decade or more. 

Keywords: Roza; basalt; flood lava volcanism; magma degassing; volcanic gases; volcanic glass chemistry; volcanic aerosols 
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asteroid impacts in geologic history (e.g., Stothers 
and Rampino, 1990; Rampino and Caldeira, 1992; 
Stothers, 1993a,b). This relationship founded an idea 
that these three events may be interrelated, where 
flood basalt episodes and mass extinctions are a 
direct consequence of impact cratering. However, the 


present general consensus, supported by age relation- 
ships at the Deccan Traps (Courtillot, 1994), is that 
impact cratering and flood basalt episodes are not 
cause and effect, but rather products of separate and 
potentially catastrophic events. This view led to the 
development of two competing hypotheses: (a) that 
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Fig. 1. (a) Map of the Roza lava flow field showing locations where the samples used in this study were collected (see text for discussion). 
Also shown is the intraflow stratigraphy at measured Roza sections (see key for details). Roza vent lineament is indicated by the heavy 
broken line. Sample locations and emplacement units are designated as in Table 1. Inset: northwestern United States, showing the location 
of the Columbia River Basalt Group (heavily stippled) and the Roza flow (lightly stippled). WA = Washington; ID = Idaho; OR = Oregeon. 
(b) Maps showing the distribution of individual emplacement units within the Roza member. The first sketch is a geological map showing 
all emplacement units. The remaining maps show the inferred original distribution of each emplacement unit in stratigraphic order from El 
to E5. Modified from Martin (1991). 








I 


T.h. Thordcirson, S. Self /Journal of Volcanology and Geothermal Research 74 ( 1996) 49- 73 5 1 

j impacts were the main cause for mass extinctions Swisher et al., 1992), whereas the development of 

r (e.g., Alvarez et al., 1980, Raup, 1986; Davis et al. , the flood basalt hypothesis has been held back by a 

j 1984; Sepkoski, 1989); or (b) that mass extinctions lack of data on the potential atmospheric H 2 S0 4 - 

weie a direct consequence of flood basalts volcanism aerosol yield by flood basalt eruptions and the mode 

| (e.g., Vogt, 1972; Courtillot et al., 1988; Stothers, of atmospheric loading. 

1993a; Courtillot, 1994). The impact hypothesis has In this paper, we address the degassing history 

gained a higher degree of acceptance in recent years, and volatile budget of a single flood basalt eruption, 

due mainly to detailed studies on the Chicxulub the one that produced the - 14.7 Ma Roza member 

impact structure and its close association with the of the Wanapum Formation in the Columbia River 

K/T boundary extinction (e.g., Hildebrand et al.. Basalt Group (CRBG), Washington, USA. This is 

1991; Sigurdsson et al., 1991; Sharpton et al., 1992; accomplished by measurements of the concentrations 



Fig. 1. (continued). 
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of S, Cl and F in the glassy products from various 
phases of the eruption. The Roza lava is favored over 
other CRBG lava flows because it is one of the best 
exposed and most distinctive flows of regional extent 
within the Columbia River basalt sequence, with 
outcrops containing the products of all the major 
eruption phases (Waters, 1961; Mackin, 1961; Bing- 
ham and Grolier, 1966; Wright et ah, 1973; Martin, 
1989, 1991). The whole-rock chemistry and the gen- 
eral stratigraphic relations of emplacement units 
within the lava field are well known (Martin, 1989, 
1991) and in terms of eruption dynamics, the Roza 
member can be taken as representing a typical 
large-volume Columbia River flood basalt event 
(Swanson et ah, 1975). Earlier, Ewart (1987) re- 
ported S -concentrations in quenched glass selvages 
from dikes (950 ppm) and pillow rinds (820 ppm) 
associated with Frenchman Springs lava flows 
(Wanapum Formation), indicating that CRBG lavas 
might be characterized by relatively high retained S 
contents. 


2. Roza characteristics and eruption history 

At present the Roza lava flow field covers ~ 
40,300 km 2 (Fig. la) and has an estimated volume of 
~ 1300 km 3 (Tolan et ah, 1989). The Roza lava is 
thought to have originated from an ~ 150-km-long 
and < 8 -km- wide NNW-trending linear vent system 
(Swanson et ah, 1975), extending from Chico Junc- 
tion, Oregon, to Winona, Washington. Outcrops of 
Roza dikes and near-vent deposits delineate the orig- 
inal Roza fissure system (Fig. la), and proximal 
outcrops show that near-vent accumulations form a 
complex association of pyroclastic deposits and lava 
flows. At Winona (RWI in Fig. la), crudely bedded 
units of lapilli-size scoria, 1-10 m thick, are found 
intermingled with thin units of vuggy pahoehoe lava 
and capped by a 0.5-2. 0-m-thick unit of densely 
welded spatter that grades upwards into a 7-15-m- 
thick fountain-fed, clastogenic lava unit (Fig. 2). 
Detailed field measurements (Thordarson, 1995) 
show that the scoria beds consist of uniform fine to 
medium-size scoria lapilli with near-horizontal bed- 
ding and are largely deprived of bomb size material. 
These textural properties indicate that the Winona 
complex represents a proximal ( < 5 km from source) 


accumulation of scoria, spatter and lava, but they are 
not vent ramparts as suggested by Bingham and 
Grolier (1966) and Swanson et ah (1975). In addi- 
tion, the small dikes ( ~ 25 cm wide) found in the 
Winona pyroclastics are surface-fed from the flow 
above, dribbling into cracks in the scoria or spatter 
deposits underneath. Scoria beds, similar to those at 
Winona, are interbedded with thin pahoehoe flows 
up to a distance of 10 km to the southeast of Winona 
(Fig. 3) and may represent the same tephra fall unit. 
Further south (between RMG and RRG in Fig. la), 
the near-vent accumulations are exposed in several 
outcrops extending over a distance of ~ 70 km and 
mostly consist of welded to rheomorphic spatter 
units. 

The Roza lavas are quartz-normative tholeiitesf 
similar to the main CRBG series (i.e., the Imnaha 
Grande Ronde-Wanapum formations), with uniform 
major- and trace-element abundances (Wright et ah, 
1973; Swanson et ah, 1979). Martin (1989, 1991) 
identified six regionally extensive chemical subtypes 
within the Roza member. This classification is based 
on subtle intraflow changes in Cr, Zr, Nb, Ti0 2 
P 2 0 5 and CaO that in general correlate with the§ 
internal stratigraphy of the Roza member on a re 
gional scale. Upward increasing Cr abundances| 
within the Roza lava succession are the most pro- 
nounced chemical change measured, suggesting thai 
later-formed lava units represent slightly less evolvec 
batches of magma. Regional correlation of these 
units based on the intraflow stratigraphy and change: 
in the lava chemistry implies that the Roza lava ca 
be grouped into at least five partly overlapping em 
placement units (Fig. la and b). On the scale of 
single outcrop, an emplacement unit consists of on 
to five pahoehoe flow units that range in total thick-i 
ness from 5 to 50 m (Thordarson, 1995). Occurrenc^f 
of multiple lava units on both outcrop and regions 
scale demonstrates the compound nature of the lav 
flow field. Roza pahoehoe sheet flow units preserv 
distinctive features such as tumuli, lava inflatio 
clefts and lava rise sutures, which are typical fo| 
inflated pahoehoe lavas (Walker, 1991; Hon et al. 
1994). Systematic measurements of the vertical dis- 
tribution of crystallinity, joint pattern and vesicula 
tion features in Roza sheet flow units show that the! 
main changes in these parameters converge at twc| 
common levels within each unit, one at nor 
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heights of -0.6 and the other at -0.1 . Most 
importantly, the upper level is associated with trans- 
formation of vertical vesicle cylinders to horizontal 
vesicle sheets. The cylinders and the sheets form 
during late stage of emplacement or very shortly 


thereafter, indicating that thick coherent surface crust 
formed on each unit during the process of emplace- 
ment (Thordarson, 1995; Self et al., 1996). Such a 
surface crust is characteristic for inflating pahoehoe 
sheet flows and grows gradually during emplacement 



(b) 



Fig. 2. Exposure of Roza near-vent accumulations at Winona (location RWI on Fig. la; see also Table 1). (a) Photograph of outcrop, 
boundaries of units are outlined by dashed lines, (b) Line drawing outlining main units. Crudely bedded scoria deposit (5) interbedded with 
thm and highly vesicular pahoehoe flow units ( P ). This sequence is capped by thick clastogenic lava ( CL ); unit is partly eroded where 
photograph was taken. Clastogenic lava has also invaded cracks (cl) that had formed in the underlying scoria bed. Scale is 1 m. 
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by accretion from below. Furthermore, the total 
thickness of the crust is a measure of the emplace- 
ment time for each unit because its growth rate is a 
function of heat lost by conduction. Making a first- 
order assumption that the Roza lava developed in a 
broadly similar fashion to similarly emplaced sheet 
flows on Kilauea, Hawaii (Hon et al., 1994), an 
emplacement time may be estimated by the expres- 
sion, t = 164. 8C~, where 1 is time in hours and C is 
crust thickness in meters. Accordingly, the measured 
surface crust thickness (6.9 and 3.4 m, respectively) 
of the sheet flow units shown in Fig. 4 indicate that 
the lower unit was emplaced in 7846 hours (—11 
months) and the upper in 1905 hours (— 2.5 months). 

The time elapsed between the emplacement of 
Roza lava units is difficult to estimate. However, a 
variety of evidence indicates that the Roza lava flow 
field was produced over a relatively short period of 
time. Even though the Roza member is composed of 
multiple flow units and several emplacement units, 
contacts between these units are often obscured or 
completely obliterated by syn-eruption welding and 
assimilation of the lava, indicating that a short time 
elapsed between the formation of these units (Swan- 
son et al., 1975; Thordarson and Self, 1993b). Com- 
mon preservation ol delicate pahoehoe surface tex- 
tures, along with near- vent occurrences of undis- 
turbed scoria layers between lava units, demonstrate 
that no significant erosion occurred between the 
emplacement of successive units. Soil or other sedi- 
ment is not found between any of the Roza lava 
units, further suggesting that a short time elapsed 
between emplacement of successive lava units. Ex- 
tensive sedimentary horizons were produced between 
the formation of major lava flow fields within the 
Columbia River plateau around the time of the Roza 
eruption (Swanson et al.. 1975). 

The Roza member is thus a compound lava flow 
field, consisting of pahoehoe lava units that were 
primarily emplaced as inflating sheet flows. We 
estimate the duration of the eruption by measuring 
systematically the thickness of the surface crust of 
Roza lava units at outcrops across the Columbia 
plateau. Using the durations obtained from the above 
equation indicates that the — 1300 km 1 Roza flow 
was formed by an eruption that lasted 5-15 years 
(Thordarson. 1995). Consequently, the calculated 
mean mass eruption rates are 7.0-22.0 X Hr kg/s. 



Fig. 3. Outcrop exposure ol Roza scoria layer (S) sandwiched 
between two Roza lava units (I ) (location RVL on Fig. la). 1 he 
scoria layer is - 40 cm thick and the lower and upper boundaries 
are indicated by white broken lines. Lower halt ot the layet is 
noil-welded and consists of line to medium size lapilli scoria, the 
upper half is moderately to densely welded presumably from the 
heat provided by the overlying lava flow. Scale is 1.0 m. 

or approximately two orders magnitude lower than 
the previously estimated eruption rates of 3.0 X 10 J 
kg/s, which was based on hypothetical fluid dy- 
namic considerations (Shaw and Swanson. 1970; 
Swanson et al., 1975). 


3. Methods 

3.1. Sample locations and descriptions 

Samples used in this study are from quenched 
glasses of Roza scoria deposits, dike and lava sel- 
vages. and glass inclusions contained in phenocrysts. 
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The samples were collected at twelve different loca- 
tions representing three of the five emplacement 
units forming the Roza lava flow field (Table 1 ; Fig. 
la). The approach employed was to collect samples 
at locations that had been noted by previous workers 
as providing a complete, or near-complete, section 


through the Roza member where the intraflow 
stratigraphy is known. Furthermore, we concentrated 
on sections sampled by Martin (1991) because his 
work provides detailed information on whole-rock 
major- and trace-element chemistry. Quenched glassy 
products from the Roza member are often remark- 


Structural Jointing Vesiculation Crystallinity 

component habits features 

meters 



Fig. 4. Section showing the main structural components of Roza lava at Banks Lake, Washington (location RBL on Fig. la) where it 
consists of two sheet flow units. The main structural components, the lava crust, lava core and basal crust, are shown on the left. The scale 
gives the thickness in meters. The left side of the stratigraphic column shows the jointing pattern of the lava. The right side of the column 
indicates vesiculation features. VZ denotes a vesicle zone, prefixes B = basal and T= top. MV and HVS indicate a horizon of megavesicles 
and horizontal vesicle sheets (equivalent to segregation veins; Goff, 1996), respectively. Vesicle cylinders ( VC ) are found as traces of 
discontinuous pipes in the lava core. Also shown are the vertical changes in crystallinity within each unit, ranging from hypohyaline ( Hhy ) 
through hypocrystalline ( Hhc ) to holocrystalline (Hoc). The boundary between the lava crust /lava core and the basal crust/lava core are 
indicated by dashed lines. Modified from Thordarson (1995). 


56 


T.h. Thordarson, S. Self / Journal of Volcanology and Geothermal Research 74 (1996) 49-73 


Table 1 


Roza sample locations and sample descriptions 


Sample location 

Sample # 

Description 

Emplace- 
ment unit 

Designation 
of Martin 
(1991) 

Dike selvages 

Roadcut, Rattlesnake Grade, Asotin County, WA 

RRG-D2-C8 

Glass selvage from E 

R-El 

R-IA 

SW1/4, SE1/4. s.25, T.7N., R.44E (45°04'N; 117°10'W) 
Elevation 635 m, dike width 1 .0 m, orientation 006°, Dip 82°E 

RRG-D2-C10 

margin of dike (C8, CIO) 
Single injection 

R-El 

R-IA 

Roadcut, Rattlesnake Grade, Asotin County, WA 

RRG-D3-C3 

Glass selvage from W and E 

R-El 

R-IA 

SW1/4, SE1/4. s.25, T.7N., R.44E (45°04'N; 117°10'W) 

RRG-D3-C5 

margin of dike (C3, C5, C6) 

R-El 

R-IA 

Elevation 635 m, dike width 2.2 m, orientation 358°, Dip 87°E 

RRG-D3-C6 

Two injection events 

R-El 

R-IA 

Roadcut, Rattlesnake Grade, Asotin County, WA 

RRG-D8-C13 

Glass selvage from W and E 

R-El + E2 

R-IA + IIA 

SW1/4, SE1 /4. s.25, T.7N., R.44E (45°04'N; 117°10'W 
Elevation 635 m, dike width 4.7m, orientation 359°, Dip 90° 

RRG-D8-C15 

margin of dike (Cl 3, Cl 5) 
Two injection events 

R-El + E2 

R-IA + IIA 

Roadcut, Casey Creek, Garfield County, WA 

SW1/4, SE1/4. s.35, T.14N., R.42E (46°1 l'N; 117°35'W) 

Elevation 480 m, dike width 2.5m, orientation 330°, Dip 85°W 

RMG-D1-C63 

Glass selvage from W 
margin of dike 
Single injection 

R-E2 

R-IIA 

Proximal scoria deposits 

Roadcut, Endicott West Road, Whitman County, WA 
SW1/4, SW1/4. s.22, T.17N., R.40E (46°57'N; 117°50'W) 
Outcrop consists of 2 Roza lava units separated by scoria 
deposits 

RWI-C60 

Glass selvage from welded 
contact of scoria deposit and 

base of the top flow unit 

R-E5 

R-IV 

Roadcut, Endicott SW Road, Whitman County, WA 
SW1/4, SE1/4. s.13, T.16N., R.40E (46°51'N; 117°45'W) 
Outcrop consists of 2 Roza lava units separated by a scoria 
layer 

RVL-C46 

Glass from welded scoria 
deposit 

R-E4? 


Lava selvages 

Summer Falls, Grant County, WA 

RSF-C26 

Glass selvage from base of 

R-El 

R-IA 

NW 1/4, NE1/4. s.15, T.23N., R.28E (47°50'N; 119°28'W) 

RSF-C41 

the lowest flow unit 

R-El 

R-IA 

Outcrop consists of 3 Roza flow units 

RSF-C37 

(C26, C41) Glass selvage 

R-E2 

R-IIA 


RSF-C40 

from base of middle flow 

R-E2 

R-IIA 


RSF-C42 

unit (C37, C40, C42) 

R-E2 

R-IIA 

Banks Lake, Grant County, WA 

SW1/4, SW1/4. s.18, T.22N., R.26E (47°72'N; 119°33'W) 
Outcrop consists of 2 Roza flow units 

RBL-C56 

Glass selvage from base 
of lower flow unit 

R-E2 

R-IIA 

Soap Lake, Grant County, WA 

SW1/4, NE1 /4. s.8, T.26N., R.28E, (47°25'N; 1 19°28'W) 
Outcrop consists of 3-5 Roza flow units 

RSL-C43 

Glass selvage from base of 
upper flow unit 

R-E2? 


Frenchman Springs Coulee, Grant County, WA 
NE1/4, NE1 /4. s.29, T.18N., R.23E (47°03'N; 119°56'W) 

RFS-C44 

Glass selvage from base 
of lower flow unit 

R-E2 

R-IIA 

Outcrop consists of 2 Roza flow units 

RFS-C53 

Glass selvage from inflation 
cleft within lower flow unit 

R-E2 

R-IIA 

Devils Canyon, Franklin County, WA 

RDC-C50 

Glass selvage of lower 
flowunit 

R-E4 

R-III 
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Table 1 (continued) 


Sample location 


SE1/4, NW1/4. s.9, T.13N., R.34E (46°63'N; 118°56'W) 
Outcrop consists of 2 Roza flow units 


Sample # Description Emplace- Designation 

ment unit of Martin 
(1991) 

RDC-C48 Glass selvage from large R-E4 R-III 

surface breakout on lower flow 
unit 

RDC-C49 Rootless cone deposit at R-E4 R-III 

base of lower flow unit 


Roadcut east of Lyle, Klickitat County, WA RLG-C55 Glass selvage from base 

NW1/4, SW1/4. s.35, T.3N., R. 12E (45°7fN; 121°28 , W) of lowest flow unit 

Outcrop consists of 2-3 Roza flow units 


Crystalline lava 

Devils Canyon, Franklin County, WA RDC266916 Crystalline lava core of 

SE1/4, NW 1/4. s.9, T.13N., R.34E (46°63'N; 118°56'W) lower flow unit 

Outcrop consists of 2 Roza flow units 


Yakima River, Yakima County, WA RDC236912 Crystalline lava core of 

NW1/4, SE1/4. s.9, T.14N., R.19E (46°43'N; 120°27'W) flow unit 

Outcrop consists of a single Roza flow unit with small lava 
lobes at its base. 


R-E2 R-IIA 


R-E4 R-III 


R-E2 R-IIA 


ably pristine considering its ~ 15Ma age, but com- 
monly are partly devitrified with patches of fresh 
glass preserved. 

Between 0.5 and 1.0 kg of rock was collected at 
each location and later examined under a binocular 
microscope to select fresh glass for thin sections and 
microprobe analyses. Four Roza feeder dikes were 
sampled (Table 1) and where possible, the glass 
selvage of both dike margins was collected. The 
composition of these dikes implies that they were 
feeders to emplacement units El and E2 (Martin, 
1989, 1991). Proximal tephra and spatter deposits 
were collected from five locations, but only two 
samples were found to contain usable glass for mi- 
croprobe analyses. Glassy basal lava selvages were 
collected from proximal to distal locations, but good 
glass for microprobe analyses was found only in the 
samples from one proximal location (sample RWI- 
C55) and six medial and distal locations (Table 1). 
The lava samples are of emplacement units El, E2 
and E4, and are representative of the spatial and 
inferred temporal distribution of the Roza lava flow 
field. Lava selvage samples from proximal localities 
are under-represented in this study because only a 
few outcrops provide exposure of the base of flow 
units in the proximal part of the Roza lava flow field. 


The Roza lava is the uppermost member of the 
porphyritic lava sequence of the Wanapum Forma- 
tion and is characterized by a distinctive phenocryst 
population, dominated by 0.3 to 1.0 cm long single, 
lath-like, honey-colored plagioclase crystals. Modal 
analysis by counting > 1000 points in several repre- 
sentative thin sections gave a mean macrophenocryst 
abundance of 9.2 + 3.0 vol.%, with plagioclase 
(An 54 _ 68 ; 6.4 vol.%), olivine (1.6 vol.%), clinopy- 
roxene (1.2 vol.%) and a trace of FeTi oxides. The 
macrophenocrysts occur as single crystals and as 
glomerocrystic aggregates consisting of a single min- 
eral phase or a mixture of the mineral types listed 
above (Fig. 5a). Where the four phases occur to- 
gether in a cluster, the order of crystallization is 
generally plagioclase, then olivine and/or clinopy- 
roxene, followed by FeTi-oxides. All of the 
macrophenocryst phases contain glass inclusions, ei- 
ther as clear light brown glass or brown to dark 
brown glass (Fig. 5b). Glass inclusions are most 
abundant in plagioclase and olivine, and are nor- 
mally of near- spherical or oblate shape with a com- 
mon size of 40-200 gm. The groundmass in crys- 
talline Roza lavas has seriate intergranular texture 
with crystal sizes ranging from <0.1 mm to 1.5 
mm. The principal groundmass mineral phases are 
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Fig. 5. (a) Transmitted-light photomicrographs of Roza dike selvage (sample C-15) demonstrating the holohyaline and vesicular nature of 
the groundmass glass. The large glomerocryst, consisting of plagioclase, clinopyroxene and FeTi oxide, is surrounded by a cloak of brown 
to dark brown non-vesicular glass. Stringers of brown glass extending from the cloak into the surrounding groundmass suggest that the 
brown glass was undergoing partial assimilation when chilling occurred. Numbers show measured sulfur concentration (ppm) of the 
groundmass glass at indicated spots. Bar is 200 yarn. (b) Transmitted-light photomicrographs of glass inclusions in phenocrysts in Roza dike 
selvage (sample C-6). Several <150 /xm brown to dark brown glass inclusions occur in plagioclase. Notice the dark brown glass stringers 
and inclusions in the groundmass glass. Bar is 200 jxm. 

plagioclase, augite, olivine, magnetite and ilmenite. 

The groundmass also contains trace amounts of ap- 
atite and interstitial brown glass. 

Dike selvages exhibit a hyalo-ophytic texture 
where macrophenocrysts ( — 10 vol.%) are embedded 
in a light brown sideromelane groundmass (Fig. 5a 
and b). The groundmass glass is moderately vesicu- 
lar, with two distinct populations, scattered 
macro vesicles (0.5-3 mm; < 3 vol.%) and a much 
greater abundance of micro vesicles ( < 0.1 mm; ~ 10 
vol.%). Overall vesicularity of the dikes ranges from 
10 to 25 vol.%, and in the dikes at Rattlesnake Grade 
(RRG in Fig. la) the vesicles are normally concen- 


trated in 1.0- to 10.0-cm-thick vertically oriented 
bands distributed evenly across the dike. These bands 
normally thicken toward the center of each dike and, 
correspondingly, the size of vesicles increases. Vesi- 
cles are usually spherical and < 1.5 cm in diameter. 

Many phenocrysts of the dike selvages are sur- 
rounded by a cloak of brown to dark brown, non- 
vesicular glass only found in the dike samples (Fig 
5a), which has a different volatile chemistry than the 
surrounding groundmass glass but the same major- 
element abundances (see below). The boundary be- 
tween the sideromelane groundmass glass and the 
brown glass is sharp but irregular because the browr 
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Fig. 6. Transmitted-light photomicrograph of Roza lava selvage (sample C-26) demonstrating the hyalo-ophitic nature of the sample. Bar is 
400 /xm. 


glass is also found as stringers being sheared out 
from these cloaks by flow related stresses (Fig. 5b). 
It also occurs as irregular- shaped inclusions in the 
groundmass glass. 

Roza tephra is composed of fine grained to 
medium size scoria lapilli, characterized by a uni- 
form grain size distribution and a moderate to low 
abundance of fines (< 1 mm fraction is < 15 wt.%). 
The scoria clasts are moderately to highly vesicular 
(35-80 vol.%) often with fluidal or fused surfaces, 
typical of achneliths (Walker and Croasdale, 1973). 
Crystal content of the scoria is typically in the range 
of 2-12 vol.%. All of the scoria samples collected 
show evidence of post-emplacement alteration (de- 
vitrification and/or palagonitization), although the 
degree to which they have been affected varies con- 
siderably and occasionally the cores of clasts contain 
some reasonably fresh glass. In outcrop, Roza scoria 
deposits are covered by later lava flow units and the 
top 20-100 cm of the scoria are commonly partly to 
densely welded, presumably due to heating from the 
overlying flow (Fig. 3). These welded scoria deposits 
often contain unaltered glass. 

The lava selvages are hypocrystalline with ~ 30 
vol.% crystals (phenocrysts and groundmass miner- 
als), wherein the seriate textured mineral assemblage 
of plagioclase, olivine and clinopyroxene (±FeTi- 
oxide) is embedded in light brown or dark brown 


sideromelane glass or black opaque tachylite glass 
(Fig. 6). The crystallinity of the lava selvages is two 
to three times higher than in quenched dike and vent 
products, suggesting that ~ 20 vol.% of the melt 
crystallized during lava emplacement (Thordarson, 
1995). 

3.2. Analytical procedure 

Modes and vesicularity of representative tephra 
and lava samples were estimated by counting ~ 1000 
points per thin section. Microprobe analyses were 
obtained at University of Hawaii using instrumental 
settings of 15 kV accelerating voltage, 15-20 nA in 
a beam-current-regulated mode and a 5-20 /a m di- 
ameter focused beam for major-element analyses of 
glass. Analyses of S and Cl were conducted by the 
Cameca Trace program, using the same conditions 
and counting for 200 s. Precision (2cr) estimated 
from counting statistics is 80 ppm for S and 60 ppm 
for Cl. The same samples where also analyzed for S, 
Cl and F by using the CSIRO-trace routine of Robin- 
son and Graham (1992) with 15 kV accelerating 
voltage, 80 nA beam current, 10 /xm focused beam 
and a counting time of 400 s. Estimated precision 
(2a) is 35 ppm for S, 30 ppm for Cl and 90 ppm for 
F. Raw data were corrected using the standard ZAF 
procedures and all standard deviations are calculated 
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as 2a. Microprobe analytical procedures are dis- 
cussed further in Thordarson et al. (1996). Whole- 
rock concentrations of S, Cl and F in two samples 
from the crystalline lava core were measured at the 
Activation Laboratories Ltd. in Ancaster, Ontario, 
Canada. 


4. Chemical analyses, results and implications 

4.1. Inclusions 

Major-element microprobe analyses of 42 glass 
inclusions in Roza phenocrysts were performed to 
locate inclusions with compositions similar to the 
host Roza lava (Table 2) and to identify glass inclu- 
sions that best represent the pre-eruption melt com- 
position. Several inclusions have a slightly more 
primitive composition, whereas others showed clear 
evidence of extensive post-entrapment crystallization 
of the host phenocryst or contained sulfide globules 
and were excluded from the study. A total of four- 
teen glass inclusions were found to exhibit major- 
element abundances (Mg# = 33.81 ± 0.81) that 
closely match the averaged bulk composition of the 
Roza member (Table 2) and these are taken as 
representative of the pre-eruption melt composition. 
The small difference found between Roza glass in- 
clusion and bulk-rock compositions is consistent with 
a pre-eruption crystal content of ~ 9 vol.% of the 
macrophenocryst assemblage described above. 

Sulfur concentrations in the selected glass inclu- 
sions fall between 1750 and 2100 ppm (mean 1965 
+ 110 ppm; Table 2) and comparable to those found 
in inclusions in phenocrysts from the 1783-1784 
A.D. Laki basaltic fissure eruption in Iceland 
(Thordarson et al., 1996). Chlorine concentrations 
vary between 200 and 350 ppm (mean 295 + 65 
ppm) and fluorine values fall between 1200 and 
1450 ppm (mean 1310+ 110 ppm). Sulfur concen- 
trations in the inclusions are in fair agreement with 
the S-FeO solubility relations in tholeiitic melts 
(Mathez, 1976; Wallace and Charmichael, 1992), 
suggesting that prior to eruption the Roza melt was 
near sulfur saturation (Fig. 7). 

Chlorine abundances in Roza inclusions fall within 
the range of 200—500 ppm, similar to values that are 
commonly reported for evolved tholeiites from other 


regions of basaltic volcanism (e.g., Sigvaldason and 
Oskarsson, 1986; Michael and Schilling, 1989, Gar- 
cia et al., 1989; Metrich, 1990). However, fluorine 
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Fig. 7. a) S (ppm) variation as a function of FeO (wt.%) in Roza 
samples. O = Roza type inclusions; O = dike selvages; A = 
scoria clasts; A = lava selvages; • = mean of crystalline lava 
(whole-rock). Solid line outlines the compositional field of mid- 
ocean ridge basalts (MORB) as defined by the data of Mathez 
(1976), Aggrey et al. (1988) and Wallace and Charmichael (1992). 
Coarse dashed line outlines the compositional field of Kilauea 
melt inclusions, using the data of Hams and Anderson (1983) and 
Anderson and Brown (1993). Fine dashed line shows composi- 
tional field of melt inclusions from basaltic flood lava eruptions in 
Iceland, using the data of Metrich et al. (1991), Thordarson (1995) 
and Thordarson et al. (1996). The composition of Roza melt 
inclusion are at the upper boundary of the MORB field, suggest- 
ing that the magma was near sulfur saturation prior to eruption 
(Wallace and Charmichael, 1992). Arrows indicate the trend of 
magma degassing during the Roza eruption, (b) S (ppm) variation 
as a function of FeO (wt.%) in Laki samples. O = Laki type 
inclusions; O = phreatomagmatic tephra; A = strombolian tephra; 
a = lava selvages; # = mean of crystalline lava (whole-rock). 
Arrows indicate the trend of magma degassing during the Laki 
eruption. The fine dashed line outlines the compositional field of 
melt inclusions from basaltic flood lava eruptions in Iceland. 
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abundances in Roza are unusually high or two to five 
times higher than reported for most evolved tholei- 
ites. The F/Cl value for the Roza magma is 4.4. 
Oceanic tholeiites appear to be characterized by 
F /Cl-values in the range of 0.5- 1.5; Pacific MORBs 
normally have F/Cl < 1 but Atlantic MORBs have 
F /Cl > 1 (Schilling et al., 1980; Michael and 
Schilling, 1989). Metrich (1990) showed that Etnean 
tholeiites are enriched in Cl and F compared to 
MORB, and characterized by F/Cl = 0.5- 1.0. 
Evolved quartz tholeiites in Iceland are typified by 
relatively high F/Cl- values of 1-3, which is at- 
tributed to assimilation of F-rich crustal material by 
magmas residing in chambers or reservoirs in the 
crust or at the crust-mantle boundary (Sigvaldason 
and Oskars son, 1986). 

Martin (1991) demonstrated that the enriched con- 
centrations of large ion-lithophile-elements in Roza 
lava (e.g., Ba = 600 ppm, Rb = 30.6 ppm, K = 1.26 
wt.%) cannot be due to crystal fractionation of man- 
tle melts and argued for minor contamination of the 
magma by “granitic” crustal material to explain 
their elemental abundances. The elevated fluorine 
and high F /Cl value of Roza glass inclusions also 
favor this interpretation and suggest assimilation of a 
Ba- and F-rich material. 

4.2. Dike selvages 

The major-element composition of groundmass 
glasses from dike selvages is very uniform (Mg# = 
33.38 ± 0.46) and comparable to that of the glass 
inclusions, with the exception that Na 2 0 concentra- 
tions are slightly lower in the inclusions (Table 2). 
The same applies to the major-element concentra- 
tions of the brown to dark brown cloaks of non- 
vesicular glass that often surround macrophenocryst 
clusters in the dike samples (Table 3; see also Fig. 
5a). It is noteworthy that the major (and volatile) 
element abundances of glass closest to the corre- 
sponding phenocryst, or within a distance of 60 /xm 
of the crystal (samples glass cloaks-tal in Table 3), 
most closely resemble those of the inclusions. Both 
glasses exhibit lower Na 2 0 concentrations than the 
groundmass glass and also lower totals (E ~ 98.10). 
These lower totals cannot be attributed to analytical 
procedure, because precisely the same results were 
obtained when analyses of the same points were 
repeated after re-polishing the sample. 


The volatile concentrations in the dike selvages 
are depleted compared to the inclusions, but they 
also demonstrate an interesting pattern. The S values 
of the light brown sideromelane groundmass glass 
range from 990 to 1205 ppm (mean 1110 + 90 ppm; 
Fig. 7; Table 3), Cl ranges from 206 to 280 ppm 
(mean 245 ± 30 ppm) and F values vary between 
990 and 1050 ppm (mean 1020 ± 25 ppm). On the 
other hand, the non-vesicular brown glass surround- 
ing the macrophenocrysts has higher S, Cl and F 
concentrations than the groundmass glass but similar 
to those in the glass inclusions (Fig. 5a; Table 3), 
especially the glass closest to the crystals (glass 
cloaks-tal in Table 3). The dark brown glass stringers 
and inclusions found in the groundmass glass have 
volatile concentrations that are intermediate between 
these two end members. Microscopic observations 
showed that the brown glass contains no sulfide- 
bearing microlites that could account for the higher 
sulfur abundances and the similarity in the major- 
element composition of the glasses rules out volatile 
enrichment due to local changes in the magma com- 
position caused by crystal growth. These cloaks are 
distinctive structures that only occur around 
macrophenocrysts in quenched Roza dike selvages 
and are in the state of being dispersed into the 
surrounding magma. They lack vesicles and are char- 
acterized by high sulfur abundances that decrease 
sharply with distance from the macrophenocryst 
clusters. These observations indicate that the glass 
cloaks represent an almost pristine melt fraction that 
was attached to the phenocrysts and was carried up 
towards the surface without exsolving much of its 
sulfur. Thus, it is possible that the low totals ob- 
tained for both the inclusions and the cloaks are from 
magmatic C0 2 and H 2 0 dissolved in the melt (Sis- 
son and Layne, 1993; Sisson and Grove, 1993), 
suggesting that prior to eruption the Roza magma 
may have contained upto 1. 5-2.0 wt.% C0 2 and 
H 2 0 in solution. In addition to glass inclusions, the 
cloaks provide further valuable information on the 
state of the Roza magma at depth prior to eruption. 

4.3. Scoria, glassy lava selvages and crystalline lava 

Major-element concentrations of the groundmass 
glass from the welded scoria are similar (Mg# = 
32.98 ± 0.42) to the glass inclusions, dike selvages 


Table 2 

Representative chemical analyses of Roza eruption products 


Sample type Si0 2 TiQ 2 A1 2 0 3 FeO MnO MgO CaO Na,0 K,0 


Whole-rock a 
Std. 

Inclusions 

Std. 

Dike selvages 
Std. 

Vent deposit 
Std. 

Lava selvages 
Std. 


51.00 

0.27 

50.52 

0.39 

51.45 

0.29 

51.43 

0.02 

51.47 

0.42 


3.17 

0.06 

3.37 

0.09 

3.40 

0.02 

3.31 

0.10 

3.86 

0.17 


13.63 

0.24 

12.54 

0.28 

12.80 

0.19 

13.05 

0.18 

12.08 

0.33 


13.93 

0.40 

14.51 

0.48 

14.46 

0.17 

14.96 

0.32 

15.73 

0.49 


0.22 

0.01 

0.24 

0.05 

0.25 

0.01 

0.21 

0.01 

0.26 

0.03 


4.55 

0.18 

4.16 

0.21 

4.07 

0.12 

4.13 

0.01 

3.51 

0.44 


8.76 

0.19 

8.16 

0.33 

8.32 

0.10 

7.81 

0.12 

7.87 

0.58 


2.78 

0.19 

2.14 

0.38 

2.73 

0.08 

1.84 

0.13 

2.21 

0.56 


1.28 

0.10 

1.46 
0.18 
1.36 
0.06 
1.44 
0.09 

1.47 
0.30 




0.68 

0.03 

0.79 

0.07 

0.75 

0.08 

0.78 

0.02 

0.91 

0.09 


Cl 


0.0195 

0.0135 

0.1965 

0.0110 

0.1110 

0.0090 

0.0665 

0.0075 

0.0640 

0.0035 


0.0100 

0 

0.0295 

0.0065 

0.0245 

0.0030 

0.0175 

0.0005 

0.0234 

0.0035 


0.0830 

0 

0.1310 

0.0110 

0.1020 

0.0025 

0.0950 

0.0020 

0.1095 

0.0055 


100.00 

0.22 

98.12 
0.86 

99.64 

0.24 

99.13 
0.27 

99.47 

0.57 


Mg# 


36.79 

1.62 

33.81 

0.81 

33.38 

0.46 

32.98 

0.42 

28.37 

3.05 




N 


289/2 

35/25 

255/83 

4 

39/50 


Table 3 


Composition of different glass types fo und in Roza dike selvages 

Sample type SiP 2 TiQ 2 A1 2 Q 3 FeO MnO MgO CaO Na 2 0 K#5 


Groundmass glass 

51.45 

3.40 

12.80 

Std. 

0.29 

0.02 

0.19 

SI glass-tc 

51.29 

3.49 

12.59 

Std. 

0.35 

0.04 

0.09 

SI glass-tb 

51.35 

3.34 

12.85 

Std. 

0.28 

0.06 

0.08 

Glass cloaks-ta 

51.51 

3.34 

12.90 

Std. 

0.05 

0.06 

0.04 

Glass cloaks-tal 

50.99 

3.29 

12.70 

Std. 

0.26 

0.13 

0.09 


14.46 

0.25 

4.07 

8.32 

2.73 

1.36 

0.17 

0.01 

0.12 

0.10 

0.08 

0.06 

14.85 

0.23 

3.83 

8.14 

2.53 

1.40 

0.19 

0.03 

0.08 

0.08 

0.25 

0.02 

14.49 

0.24 

4.03 

8.24 

2.62 

1.37 

0.19 

0.03 

0.10 

0.16 

0.10 

0.05 

14.46 

0.26 

3.97 

8.27 

2.49 

1.42 

0.21 

0.04 

0.25 

0.19 

0.14 

0.08 

14.23 

0.24 

4.06 

8.28 

2.24 

1.35 

0.28 

0.03 

0.08 

0.15 

0.21 

0.05 


P2O5 S Cl F £ 


0.75 

0.08 

0.81 

0.05 

0.74 

0.06 

0.71 

0.09 

0.70 

0.12 


0.1110 

0.0245 

0.1020 

99.64 

0.0090 

0.0030 

0.0025 

0.24 

0.1255 

0.0250 

0.1110 

99.16 

0.0070 

0.0035 

0.0080 

0.79 

0.1380 

0.0300 

0.1040 

99.54 

0.0125 

0.0055 

0.0125 

0.55 

0.1795 

0.0240 

0.1115 

99.80 

0.0035 

0.0020 

0.0140 

0.39 

0.1881 

0.0280 

0.1170 

98.20 

0.0090 

0.0040 

0.0140 

0.60 


Mg# 

N 

33.38 

0.46 

255/83 

31.47 

0.46 

8 

33.15 

0.36 

5 

33.25 

0.46 

39/15 

33.68 

0.60 

12/4 


-t/ represel^dark bro^gll^ th^licurs^s^tringem and ^rT maj ° r . /V ° Iatile element anal >' ses; Std - = standard deviation (2a). SI glass-tb and 

glomerocrysts; ta = analyses at distances between 60 and 200 from crystals” ta”- ^“ysesTt LTanc” < nspres ' m dark brown - "°"-vesicular glass surrounding 
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and bulk lava, whereas the lava selvages are more 
evolved (Mg# = 28.37 ± 0.35) and have somewhat 
more variable composition (Table 2). Compared to 
bulk-rock compositions, the Roza lava selvages are 
enriched in Ti0 2 (0.5 wt.%) and FeO (1.0-1. 5 wt.%). 
This differentiation of the groundmass glass is at- 
tributed to in situ crystallization (or equilibrium crys- 
tallization at a constant total composition) of the lava 
during emplacement. Roza lava selvage composi- 
tions show no correlation between the degree of 
differentiation and distance from source. 

The groundmass glass of a few lava selvages and 
scoria clasts contain a little lower CaO (~ 7.5 wt.% 
versus 8.32 wt.%) and Na 2 0 (~ 1.85 wt.% versus 
2.73 wt.%) than the dike selvages, although similar 
elemental abundances are found in the glass inclu- 
sions. The lower CaO and Na 2 0 content is not likely 
to be caused by palagonitization because K 2 0, P 2 0 5 , 
Ti0 2 and FeO, which have similar mobility to CaO 
and Na 2 0 during palagonitization (Fumes, 1978), 
are unaffected. However, the lava and scoria samples 
with low CaO and Na 2 0 glasses have more plagio- 
clase crystals (phenocrysts + groundmass mode = ~ 
22 vol.%) than the high CaO and Na 2 0 samples 
(modal plagioclase closer to ~ 15 vol.%). Thus, the 
lower abundances of CaO and Na 2 0 may have 
resulted from crystallization of more plagioclase in 
some lavas than in others. 

Sulfur values in the scoria groundmass glass have 
a mean of 665 ±75 ppm, which is comparable to the 
S content of the lava selvages (640 ± 65 ppm), but 
well below that of the dike selvages (Fig. 7; Table 
2). Chlorine concentration averages 175 + 5 ppm in 
the scoria but in the lava selvages the mean value is 
slightly higher at 245 ± 50 ppm. Similarly, the F- 
values are 845 ±150 ppm for the scoria and 1085 ± 
85 ppm in the lava selvages. Because ~ 10-20% 
j crystallization occurred during lava emplacement, 

! the concentrations of S, Cl and F are enriched com- 
pared to the original fraction of magma erupted 
| (Thordarson et al., 1996). These elements were in- 
compatible with crystallizing phases in the Roza 
magma, therefore it is expected that S, Cl and F 
abundances gradually increased in the liquid lava 
during emplacement. This volatile enrichment pro- 
cess was complicated by simultaneous loss of these 
elements through degassing (especially S) and by the 
increased solubility of S in the melt due to FeO 


Table 4 

Corrected S, Cl and F values for Roza lava selvages 


Lava selvage 
(sample) 

S 

S a 

Cl 

Cl a 

F 

F a 

RSF-C26 

0.0625 

0.0500 

0.0260 

0.0210 

0.1015 

0.0810 

RSF-C41 

0.0705 

0.0540 

0.0290 

0.0225 

0.1150 

0.0885 

RSF-C37 

0.0655 

0.0530 

0.0320 

0.0260 

0.1145 

0.0925 

RSF-C40 

0.0580 

0.0485 

0.0205 

0.0170 

0.1065 

0.0890 

RSF-C42 

0.0655 

0.0525 

0.0208 

0.0165 

0.1095 

0.0875 

RBL-C56 

0.0595 

0.0505 

0.0225 

0.0190 

0.1120 

0.0950 

RSL-C43 

0.0630 

0.0535 

0.0255 

0.0215 

0.1120 

0.0950 

RFS-C53 

0.0670 

0.0560 

0.0210 

0.0175 

0.1145 

0.0955 

RFS-C44 

0.0675 

0.0525 

0.0225 

0.0175 

0.1160 

0.0905 

RDC-C50 

0.0630 

0.0515 

0.0212 

0.0175 

0.0990 

0.0810 

RDC-C48 

0.0645 

0.0525 

0.0192 

0.0150 

0.1115 

0.0905 

RDC-C49 

0.0600 

0.0460 

0.0220 

0.0170 

0.1025 

0.0790 

RLG-C55 

0.0660 

0.0575 

0.0220 

0.0190 

0.1090 

0.0945 

Average 

0.0640 

0.0520 

0.0234 

0.0190 

0.1095 

0.0890 

Std. 

0.0035 

0.0030 

0.0035 

0.0030 

0.0055 

0.0055 


Std. = standard deviation (2a). 
a Corrected values. 


enrichment. At any rate, to find the concentration of 
these elements in the bulk lava (the original amount 
of liquid lava), it is necessary to correct for the 
amount of crystallization that occurred during em- 
placement. The corrected volatile values for the lava 
selvages are 520 ± 30 ppm S, 190 ± 35 ppm Cl and 
890 ± 55 ppm F and still show a remarkable uni- 
formity (Table 4). 

The volatile concentrations of the scoria and the 
lava selvages are very similar, and the differences 
are of the same order as the analytical uncertainties. 
The uniform volatile concentrations of the lava sel- 
vages, however, suggest that minimal melt degassing 
occurred during lava transport and that the Roza 
lavas contained ~ 520 ppm S when they flowed 
away from the vents. Using the results of Haughton 
et al. (1974) and eruption temperatures of 1140°C 
(Schiffman and Lofgren, 1982), the estimated sulfur 
capacity of the Roza lava as it flowed away from the 
vents is ~ 350 ppm. This value is lower but compa- 
rable to the obtained bulk sulfur concentrations in 
the groundmass glass of Roza scoria and lava sel- 
vages, suggesting that the vent degassing of S ap- 
proach an equilibrium stage. Finally, the S, Cl and F 
concentrations in crystalline lava are 195 ppm S, 100 
ppm Cl and 830 ppm F (Table 2). 
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5. Magmatic volatiles and degassing history 

5.1. Model for estimating volatile release by flood 
lava eruption 

Several studies (e.g., Devine et al., 1984; Andres 
et al., 1989; Thordarson et al., 1996) have clearly 
demonstrated the validity of a petrologic approach to 
evaluate the potential atmospheric sulfur yield from 
basaltic fissure eruptions. However, the applicability 
of the method may depend on the volcano type and 
pre-eruption conditions in the magma source region. 
The petrologic method appears to underestimate the 
S0 2 flux from some eruptions at convergent-plate 
margins because (a) excess sulfur is derived from 
unerupted magma or intrusions at shallow depths 
within the volcano (Rose et al., 1982: Andres et al., 
1991), and (b) gas-saturated melts residing in crustal 
magma chambers produce a separate vapor phase 
which carries the bulk of the sulfur during an erup- 
tion (Wallace and Gerlach, 1994; Gerlach et al., 
1994). These conditions do not apply to the fissure- 
fed Roza flood lava eruption because of the mafic 
nature, the large volume and the homogenous com- 
position of the products which imply that the magma 
was erupted from a compositionally uniform reser- 
voir located at the base of the crust (Hooper, 1984; 
Catchings and Mooney, 1988). Furthermore, the glass 
inclusions have similar major-element abundances as 
the degassed eruption products, but their volatile 
concentrations are much higher (Table 2), suggesting 
that they are representative of the pre-eruption Roza 
magma. In addition, a similar study of the Laki 
1783-1984 tholeiite basaltic fissure event in Iceland 
(Thordarson et al., 1996), an identical but smaller 
flood basalt eruption (Thordarson, 1995), the petro- 
logic estimates of atmospheric sulfur mass loading 
gave comparable results to independent assessments 
using information on atmospheric turbidity over Eu- 
rope (Stothers, 1996). Therefore we conclude that 
the petrologic method gives representative estimates 
for the S, Cl and F degassing of the Roza magma 
upon eruption. In applying the petrologic method, 
however, we assume that the volatile mass released 
the Roza eruption is derived entirely from the vol- 
ume of erupted magma. It is likely that unerupted 
magma contributed to the volatile mass, hence, our 
petrologic estimates should be regarded as minimum 


Eruption column and distal haze 



Fig. 8. Illustration outlining a volatile budget model for basaltic 
fissure eruptions that can be described as follows: ( 1 ) Hypotheti- 
cal magma reservoir containing magma of the mass M r (taken 
here to be equal to the erupted volume of magma) and total 
amount of dissolved volatiles m r \ the measured concentration of 
volatiles in glass inclusions, v- t , represents the concentration of 
volatiles dissolved in the magma. The total amount of volatiles in 
the magma reservoir is given by the lowermost equation. (2) Dike 
selvages containing concentration of volatiles u d . ( 3 ) Vent accu- 
mulations of pyroclastic material and lava selvages containing 
concentrations of volatiles v n which represents volatile fraction of 
m r that did not escape the magma as it emerged from the vents. 
(4) Eruption column (lava fountain plus convecting plume) con- 
taining total amount of volatiles released at the vents, m v , as 
given by the top equation. (5) Crystalline lava containing volatile 
concentration v c . (6) Amount of volatiles, m c , released by the 
outgassing of lava during and after emplacement and is given by 
the middle equation. Modified from Thordarson et al. (1996). 


values for the mass of S, Cl and F released by the 
Roza eruption. 

If the total mass of eruptives equals the amount of 
magma that effectively contributes to the volatile 
emission during an eruption, then the overall volatile 
budget of the eruption can be represented by the 
schematic model illustrated in Fig. 8. Magma de- 
gassing in effusive eruptions can be viewed as occur- 
ring in two separate stages, first at vents through 
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Strombolian- or Hawaiian-type fire fountaining and, 
later, during and after lava emplacement (Oskarsson 
et al., 1984; Thordarson et al., 1996). 

The total mass ( m xx ) of a volatile element x 
dissolved in magma residing in a reservoir at depth 
and effectively involved in the eruption can be eval- 
uated by: 
m x = M x v xi 

where M x is mass of magma erupted and u X ] is the 
mass fraction of element x in inclusions (i). The 
total mass (m x v ) of volatile element ^ that escapes 
from the magma at the vents and is dispersed into 
the atmosphere by the eruption column is estimated 
by: 

m x,y = ( V x3 ~ V x,t) M r 

where v xt is the averaged mass fraction of element 
x as measured in samples from vent deposits and 


lava selvages. Similarly, the mass released into the 
atmosphere by gaseous emissions from the lava dur- 
ing cooling and crystallization (m xc ) is estimated 
according to: 

m c = ( V x,t- V x, C ) M l 

where v xc are the mass fraction of element x in 
glassy lava selvages and crystalline lava, respec- 
tively. M y is total mass of lava. 

5.2. Application 

The amount of each volatile species carried by a 
melt to the surface is calculated by: 

ra r = M x v xi — V x pe v x i 

where V x is total volume erupted, taken here as 1300 
km 3 (Tolan et al., 1989), p is the magma density 
taken to be 2700 kg m~ 3 and v x[ is the mass fraction 


Table 5 

(a) Averaged concentration (in ppm) of volatile species in the Roza eruption products 

S Cl 

Inclusions (u;) 1965 295 

±110 ±65 

Dike selvages (y d ) 1110 245 

±90 ±30 

Scoria and lava selavages (v t ) 595 185 

±75 ±30 

Crystalline lava (v c ) 195 100 

±135 ±0 

(b) Volatile abundances (in ppm) used in mass balance calculations 


AS AC1 

Total amount released (v l - u c ) 1770 195 

Amount released in feeder dike (uj — n d ) a 855 50 

Amount released at vents (ii; — u t ) 1370 1 10 

Amount released by lava (n t — v c ) 400 85 


(c) Estimates on degree of magma degassing at various eruption stages 


S 



(%) 

(%) 

(%) 

(%) 

Total degassing (v t — v c ) 

90 

66 

37 

68 

Amount retained in solidified products ( u c ) 

10 

34 

63 

32 

Degassing in feeder dike (u ; — u d ) 

44 

17 

15 

33 

Degassing at vent (iij — t> ( ) 

70 

37 

30 

52 

Degassing by lava (u, — u c ) 

20 

29 

7 

16 


a Shows fraction of volatiles already exsolved from the magma at ~ 1 km depth in the conduit. This fraction contributes to the total amount 
of volatiles released at the vents. 
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of element x in glass inclusions, and e is the factor 
required for converting the pure element to the as- 
sumed elemental compound present in the magma 
(e.g., S0 2 , HC1 and HF). The elemental concentra- 
tions used here to calculate the amount of volatiles 
released by various stages of the Roza eruption are 
presented in Table 5. As an example, the average 
inclusion value of sulfur, 1965 ppm (0.1965 wt.%), 
transforms into ~ 14 X 10 3 megatons (Mt) of S0 2 
carried from depth to the surface by the erupting 
magma (Table 6; Fig. 9). 

The total mass of volatile element (m x h ) released 
by the eruption is taken as the difference between the 
inclusion value (m x i ) and the mass retained in the 
solidified eruption products ( m x c ). This difference 
suggests that ~ 90% of the S0 2 escaped from the 
erupted magma during the eruption (Table 5) and 
that the total mass released was ~ 12.5 X 10 3 Mt 


Table 6 

Estimates on the mass (in megatons) of volatiles dissolved in the 
Roza magma prior to eruption and released by various phases 
during the eruption 



S0 2 

HC1 

HF 

£ 

Original mass (m r ) 

13,795 

1065 

4840 

19,700 

Total mass released (m b ) 

12,425 

704 

1774 

14,903 

Mass retained in solidified 

1369 

361 

3068 

4788 

products ( m c ) 

Mass released at vents (m r — m v ) 

9617 

397 

1442 

11,456 

Mass released by lava (ra t — m c ) 

2808 

307 

333 

3448 


Total volume of magma used in the calculations of 1300 km 3 is 
from Tolan et al. (1989) and magma density is taken as 2700 
kgm -3 . 


S0 2 (Table 6; Fig. 9). The fraction of HC1 and HF 
released is considerably smaller, estimated at ~ 0.7 
X 10 3 Mt HC1 and ~ 1.8 X 10 3 Mt HF. 

The S, Cl and F concentration in the groundmass 


km 

-18 



Fig. 9. Schematic illustration of the Roza eruption (not to scale) displaying the amount of S retained in samples from various eruption stages 
and the total amount of S0 2 dissolved in the Roza magma prior to eruption. Also shown is the estimated S0 2 -yield by various eruption 
stages. The gas columns and plumes generated by the magma degassing are shaded light gray areas. Distal haze denotes the gas released at 
the vents and carried to high altitudes by the eruption columns; local haze denotes the gases rising of the lava flow. Glass inclusions 
(S = 1965 ppm ) in phenocrysts represent undegassed melt residing in the reservoir prior to eruption. Groundmass glass of dike selvages 
( mo ppm ) represent partly degassed melt that froze in the feeder dike. Groundmass glass of vent deposits ( 665 ppm) and lava selvages 
( 520 ppm) characterize melt that has gone through the stage of vent degassing, and in addition, has released fraction of its remaining 
volatiles as the lava flowed away from the fissures. Crystalline lava (795 ppm) denotes melt that has gone through the stages of vent and 
lava emplacement degassing and has lost additional fraction of its volatiles during cooling and crystallization. It also contains the residue 
volatiles that did not exsolve from the lava. 
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glass of dike selvages fall between those of the 
inclusions and vent deposits, which is consistent with 
the groundmass glass representing melt undergoing 
degassing as it was rising towards the surface. The 
elevation difference between the dike outcrops (Ta- 
ble 1) and the base of the Roza flow in nearby cliff 
exposures around Rattlesnake Grade is 600-1000 m 
(Swanson et al., 1980; S.P. Reidel, unpublished geo- 
logical map of Asotin County, 1994), indicating that 
the dike samples correspond to magma at crustal 
depths of 0.5— 1.0 km. The dike selvages represent 
rapidly chilled magma from the earliest stages of 
vertical magma injection, whereas the crystalline 
interior of the dike may represent batches from later 
stages of magma injection. Thus, the dike selvages 
provide information on the amount of volatiles ex- 
solved from the rising magma at the given crustal 
level. Roza dike selvages (Table 5) show that magma 
had already exsolved 44% of its original S and 
~ 15% of Cl and F into vesicles and/or a separate 
volatile phase rising ahead of the magma column as 
it passed crustal depths of 600-1000 m, i.e., the 
level of the Rattlesnake Grade dike outcrop. It is 
likely that the fractions of exsolved H 2 0 and C0 2 
gases were similar to that of S0 2 (Thordarson et al., 
1996). Thus, with ~ 45% of the original magmatic 
volatiles already exsolved at crustal depths of 0.5- 1.0 
km, theoretical relationships (Wilson and Head, 
1981) suggest that the Roza magma contained 1.5- 
2.0 wt.% of dissolved volatiles prior to eruption. 
This amount agrees with the estimated vapor (mag- 
matic gas) amounts based on low totals in Roza glass 
inclusions (see above), and is similar to the calcu- 
lated contents of Laki total volatiles (Thordarson et 
al., 1996). 

The S, Cl and F concentrations in the glassy 
groundmass of the scoria and lava selvages samples 
(Table 2) are taken to represent the volatile abun- 
dance in the magma as it flowed away from the 
vents. Consequently, the difference between the in- 
clusion values and the average values for the scoria 
and lava selvages is the volatile fraction released by 
the melt into the atmosphere above the vents, signi- 
fying that 70% of the total sulfur was released as the 
magma emerged from the fissures (Table 5). Thus, 
the Roza magma may have discharged ~ 9600 Mt 
of S0 2 into the atmosphere above the fissures (Table 
6; Fig. 9). The vent degassing of HC1 and HF is 


considerably less, 37 and 30%, respectively. Conse- 
quently, the eruption released about 400 Mt HC1 and 
1400 Mt HF into the atmosphere above the vents. 
The partly degassed magma flowed from the vents to 
form lava flows, either effusively and/or via foun- 
tain collapse forming clastogenic lava flows. 

The volatile-element values of the scoria and lava 
selvages (Table 4) are taken as the fraction of 
volatiles contained in the lava upon emplacement in 
the medial and distal areas of the lava field. Conse- 
quently, the concentration difference between 
volatile-element concentrations in lava selvages and 
crystalline lava denotes the amount released during 
cooling and crystallization of the lava after the flow 
had stopped its advance. Accordingly, the Roza flow 
degassed ~2800MtS0 2 , ~ 305MtHCland -330 
Mt HF during cooling and crystallization (Table 6; 
Fig. 9) or approximately one-fifth of the total volatile 
mass released (Table 5). 


6. Discussion 

6.1. Roza eruption dynamics and magma degassing 

The Roza lava field, along with other CRB Wana- 
pum lava formations, are a typical example of fis- 
sure-fed flood basalt event (Swanson et al., 1975) 
and in essence were large scale analogs of present 
day effusive basaltic volcanism. Although the vol- 
ume-ratio of pyroclastic ejecta to lava produced by 
the Roza eruption is small ( 1 /20), previous field 

studies (e.g., Bingham and Grolier, 1966; Swanson 
et al., 1975; Martin, 1991) clearly show that the 
eruption produced large amounts of tephra and spat- 
ter deposits (possibly in the range of 1-10 km 3 ). 
Outcrops show that near-vent accumulations form a 
complex association of pyroclastic deposits and lava 
flows, indicating that the Roza eruption was not only 
characterized by tremendous outpourings of lava, but 
also involved episodes of intense lava fountaining 
events with widespread dispersal of scoria. However, 
the duration of lava fountaining events during the 
Roza eruption are not known. The sulfur content of 
the Roza glasses defines a distinctive trend (Fig. 7) 
which reflects the progress of melt degassing during 
a basaltic flood lava eruption and is consistent with 
the estimate that Roza melt released — 70% of its 
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S0 2 through conduit and vent degassing processes. 
This trend is very similar to the one found for the 
Laki flood lava eruption in Iceland (Fig. 7b). With 
initial volatile content of ~ 1.7 wt.% the Laki magma 
released > 70% of the S0 2 , CO, and H 2 0 through 
vent degassing (Thordarson et al., 1996). The inclu- 
sion data suggest that the Roza magma may have 
initially contained ~ 1.9 wt.% of C0 2 and H 2 0 and 
by comparison to Laki we infer that the fraction of 
these volatile species released at the Roza vents was 
similar to that of S0 2 (i.e., ~ 70% of the total mass 
released). The large fraction of volatiles released by 
the Roza magma up on venting, along with the 
relatively small volume of pyroclastic material pro- 
duced, implies rapid and effective separation of gases 
and melt during ascent and probable development of 
a separated two-phase flow in the feeder dike. Thus, 
the Roza eruption may have been capable of main- 
taining sustained eruption columns that reached con- 
siderable altitudes in the atmosphere. 

To accurately assess the heights reached by the 
Roza eruption we need more detailed information 
than is available at present but we can make a 
general estimate. Such an estimate is heavily depen- 
dent on magma discharge and the geometry of the 
fissure system active at any given time during the 
eruption. The average mass eruption rates for Roza is 
7.0 X 10 6 -22.0 X 10 6 kg/s, but magma discharge 
during periods of peak activity in basaltic fissure 
eruptions may be six to ten times greater than the 
average (Thordarson and Self, 1993a; Thordarson et 
al., 1996). Martin (1991) suggested that individual 
emplacement units within the Roza flow field were 
produced by linear vent segments that range from 
~ 7 to ~ 75 km in length. Although magma out- 
pourings in basaltic fissure eruptions may occur along 
the entire fissure length at the beginning of an 
eruption, such outpourings are rarely sustained for 
any length of time (e.g., Wilson and Flead, 1981). 
Magma discharge is quickly localized to one or more 
much shorter segments. Indirect evidence for local- 
ization of vent activity on a Columbia River basalt 
fissure system is provided by Reidel and Tolan 
(1992), who describe a vent complex consisting of 
proximal shelly pahoehoe lavas, pyroclastic deposits 
and a lava lake typical of Hawaiian volcanism. 
Therefore we believe that the Roza eruption was 
characterized by activity similar to historic fissure 


eruptions, such as at Kilauea and Laki, except for the 
much greater magnitude of the eruption. 

If we use the mass eruption rates given above and 
assume that activity occurred simultaneously along 
the longest Roza vent segment (~ 75 km), then the 
magma discharge per meter length of fissure — 1.5 
X 10 2 kg/s. Using the model calculations of Woods 
(1993), the corresponding height of the eruption 
columns in a moist atmosphere above the lava foun- 
tains is ~ 3-6 km. Similarly, if we use the shortest 
Roza vent segment ( ~ 7 km), then the discharge per 
meter length of fissure — 2 X 10 3 kg/s and the 
eruption would have maintained 6-8.5 km-high 
eruption columns above the lava fountains. Assum- 
ing that the magma was emitted from a small fissure 
segment, taken here as a point source (discharge = 
1.7 X 10 7 kg/s), the model predicts 11-13 km for 
the height of the eruption columns. Direct observa- 
tions and model calculations indicate that the 1783 
Laki event in Iceland maintained eruption columns 
extending to altitudes of 7-13 km during the first 3 
months of activity, when peak magma discharge was 
~ 1 X 10 7 kg/s (Thordarson and Self, 1993a). A 
comparison indicates that the Roza eruption may 
have been capable of maintaining eruption columns 
that reached to similar altitudes for the estimated 
duration of the eruption (5-15 years). If the relative 
difference between the mean and the peak discharge 
during the Roza eruption were similar to that of 
Laki, then the maximum eruption rates would have 
been on the order of 10 8 kg/s. In this case the 
height of the eruption columns predicted by the 
model for the above three cases would be doubled, 
i.e., - 9, ~ 17 and ~ 28 km, respectively. This 
general treatment of the Roza eruption columns 
demonstrates that the Roza eruption was capable of 
maintaining columns that reached well into the upper 
troposphere for the entire duration of the eruption 
and had the potential to generate columns that would 
have penetrated into the stratosphere. The eruption 
column heights would have been highly dependent 
on the spatial and temporal variations in activity 
along the eruptive fissure. 

6.2. Atmospheric loading of S0 2 and possible cli- 
matic effects of Roza 

Evaluations of the potential climatic effects of 
flood basalt eruptions are strongly dependent on the 
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overall dispersal of aerosols. Because atmospheric 
transport and residence time of aerosol particles are 
different at various altitudes in the atmosphere, it is 
important to determine the height of the associated 
eruption columns from which the gases would be 
dispersed. However, we have shown that the erup- 
tion history of flood basalt events, including their 
spatial and temporal changes, can only be con- 
strained very generally at present. 

Because of the uncertainties regarding the progress 
of the Roza eruption discussed above, the following 
discussion on the possible atmospheric effects of the 
Roza eruption can only be treated in a generalized 
and somewhat speculative manner. The total amount 
of S0 2 released into the atmosphere above the Roza 
vents is estimated to have been ~ 9600 Mt (Table 6) 
and because the bulk of the degassing occurred at the 
vents much of this mass may have been lofted by the 
eruption columns to altitudes > 10 km. Plate tec- 
tonic reconstructions suggest that the paleolatitude of 
the Roza eruption site was ~ 50°N (Smith and 
Briden, 1977). If the overall structure of the atmo- 
sphere in mid-Miocene times (~ 14.5 Ma) was simi- 
lar to that of today, then the Roza columns would 
have reached the upper troposphere and lower strato- 
sphere at that latitude. 

We here make a first order attempt at estimating 
the atmospheric loading from a long-lasting flood 
basalt eruption. Assuming the eruption lasted for 10 
years, the S0 2 emissions from the Roza vents would 
have released 900 Mt per year, an order of magni- 
tude greater than the present annual anthropogenic 
S0 2 emissions on Earth (Park, 1987). This further 
implies that the annual H 2 S0 4 -yield from the S0 2 
produced by the Roza eruption was ~ 1800 Mt. This 
estimate is obtained by converting the total mass of 
S0 2 released by the eruption into H 2 S0 4 (18,000 
Mt); the estimated duration of the eruption (10 years) 
is then divided into the potential total H 2 S0 4 -yield. 
The residence time of an average aerosol particle 
(radius = 0.3 yam) in the upper troposphere and the 
lower stratosphere is on the order of 50 to 400 days 
(Jaenicke, 1984), although the atmospheric removal 
rates are also strongly dependent on the rate of 
conversion of S0 2 to H 2 S0 4 , the overall particle 
size distribution and mass of the aerosols (Pinto et 
al., 1989). Although such huge amounts of aerosols 
will fall out fast, they could be resupplied by contin- 


uous S0 2 emissions. The removal rates given above 
suggest that the Roza eruption was able to sustain 
very high atmospheric concentrations of H 2 S0 4 
aerosols, as it can be presumed that the aerosol cloud 
was being replenished at fairly steady rates (Thordar- 
son and Self, 1996). We thus estimate that Roza 
maintained an annual atmospheric H 2 S0 4 loading in 
the range of 500-1000 Mt. If the aerosol cloud was 
confined to the Northern Hemisphere, then the opti- 
cal thickness of the cloud can be estimated from: 

9 QM d 

D 16 irR 2 rp 

where r D is optical depth, M D is the aerosol mass 
(kg), R is the radius of the Earth, Q is an efficiency 
factor for scattering and absorption by the aerosol 
particles, and r and p are the aerosol particle radius 
and density, respectively (Stothers, 1984). Adopting 
typical modal values of <2 = 2, r = 3 X 10“ 7 m and 
P— 1500 kgm -3 , then the range of aerosol optical 
depths caused by the Roza aerosol cloud may have 
been 7-13. These estimates are based on current 
knowledge of aerosol formation after sudden vol- 
canic events, because we do not have models that 
account for continuous S0 2 release from sustained 
volcanic eruptions lasting for years. 

On the other hand, if we assume a global disper- 
sal of the aerosol cloud, then r D « 3—6, which is of 
similar magnitude as estimated for the 73.5 kyr Toba 
super eruption (Rampino and Self, 1993). Atmo- 
spheric perturbations of this magnitude are unprece- 
dented in recent history and the consequences of 
such massive atmospheric loading of volcanic 
aerosols are unknown and debatable, but could possi- 
bly be evaluated by recent generations of global 
circulation models. Based on the above considera- 
tions, the opacity of the Roza aerosol cloud may 
have been similar to that of an overcast sky (Fig. 
10), conditions which the eruption could have main- 
tained for ~ 10 years. Such a scenario is similar to 
that predicted for severe “nuclear winters” (Turco 
et al., 1983, 1984; Ackerman, 1988) or “volcanic 
winter” (Rampino et al., 1988; Rampino and Self, 
1993), which may be associated with immediate 
surface cooling on the order of 5-15°C. If the Roza 
eruption maintained such an aerosol cloud for a 
period of 10 years, atmospheric conditions would 
severely affect the global climate. 
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Aerosol loadings (Mt) 

Fig. 10. Fraction of sunlight transmitted through aerosol and/or fine-ash clouds of different masses (theoretical line after Turco et al., 
1984). Points refer to designated historic and prehistoric eruptions. The point for Roza refers to estimated annual mass loading of aerosols 
but the eruption may have maintained this mass loading for several years. Modified from Rampino et al. (1988). 


Indeed the Miocene paleotemperature records, re- 
constructed from oxygen and carbon isotopic compo- 
sitions in benthic and planktonic foraminifera from 
deep-sea sediments in the Pacific and Atlantic 
(Shackleton and Kennett, 1975; Vincent and Killing- 
ley, 1985; Kennett, 1986; Miller et al., 1987) and 
from the Monterey Formation in California (Flower 
and Kennett, 1993), show initiation of a distinct 
global cooling event of several degrees Celsius at 
~ 15.0-14.5 Ma or the time of the Wanapum basalts 
and Roza eruption. This cooling event lasted until 
~ 14.1 Ma. Moreover, the sharp cessation of the 
Early Miocene warming around 16 Ma that is be- 
lieved to characterize the transition from a relatively 
warm global climate to the Neogene “icehouse 
world” coincided with the peak of the Columbia 
River flood basalt activity at Grande Ronde times 
(Reidel et al., 1989). At that time the average erup- 
tion frequency and volumetric outpourings of magma 
were equivalent to one Roza eruption per 7500 years 
(Tolan et al., 1989). Although the residence time 
(years to tens of years) of the aerosol cloud from a 
single flood basalt eruption is by far too short to 
have been a direct cause for cooling events of this 
magnitude and duration, is it possible that the mas- 
sive atmospheric aerosol loading from the series of 
eruptions that comprise the Wanapum Formation, 


including Roza, served to trigger a global cooling at 
14.5 Ma? Is it possible that the whole CRBG flood 
basalt episode triggered global deterioration of cli- 
mate that eventually led to the Neogene “icehouse 
world”? 


7. Summary 

Glass inclusions in phenocrysts in the 1300 km 3 
Roza lava flow contain 1965 ppm S, 295 ppm Cl and 
1310 ppm F, and the Roza magma was near sulfur 
saturation prior to eruption. Analyses of groundmass 
glass from Roza dike selvages reveal considerably 
lower concentrations: 1110 ppm S, 245 ppm Cl and 
1020 ppm F. Groundmass glasses of scoria clasts and 
lava selvages contain 595 ppm S, 185 ppm Cl and 
920 + 75 ppm F and crystalline lava contains 195 
ppm S, 100 ppm Cl and 830 ppm F. These samples 
define a trend that represents the progress of de- 
gassing through the eruption and can be used to 
estimate the total amount of volatiles released by the j 
melt, as well as to evaluate the amount liberated by i 
various eruption phases. 

The total amount of volatiles released into the 
atmosphere by the Roza melt is estimated to be 
12,500 Mt S0 2 , 700 Mt HC1 and 1750 Mt HF. 
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Furthermore, we estimate that 9600 Mt S0 2 , 400 Mt 
HC1 and 1450 Mt HF were liberated at the vents and 
lofted by the eruption columns to 7-13 km altitudes. 
The lava is estimated to have released an additional 
2800 Mt S0 2 , 305 Mt HC1 and 330 Mt HF. These 
results indicate that the atmospheric perturbation by 
the Roza eruption may have been of the magnitude 
predicted for a severe “nuclear” or “volcanic” 
winter, with sudden and short-lived surface cooling 
of ~ 5-15°C lasting up to a decade or more. 
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